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Abstract 

An analysis of the random lattice gas in the annealed limit is presented. The statistical mechanics of 
disordered lattice systems is briefly reviewed. For the case of the lattice gas with an arbitrary uniform 
interaction potential and random short-range interactions the annealed limit is discussed in detail. By 
identifying and extracting an entropy of mixing term, a correct physical expression for the pressure 
is explicitly given. The one-dimensional lattice gas with uniform long-range interactions and random 
short-range interactions satisfying a bimodal annealed probability distribution is discussed. The model 
is exactly solved and is shown to present interesting behavior in the presence of competition between 
interactions, such as the presence of three phase transitions with different critical temperatures and the 
occurrence of triple and quadruple points. 

The lattice gas model fjj was introduced by Lee and Yang in 1952 as an application of their theory of 
condensation, due to its equivalence to the Ising model. Despite its somewhat artificial character, it has 
been extensively and successfully applied to the study of both solid-liquid and liquid-gas transitions || . It 
has also been applied to problems of adsorption of a gas on a crystal surface d|, | . 

The mapping of the two-dimensional lattice gas with nearest-neighbor interaction into the Ising model 
in the square lattice allowed the exact calculation of the coexistence curve of a fluid transition exhibited by 
the lattice gas at sufficiently low temperatures. However, the corresponding equation of state could not be 
exactly determined, due to the absence of an exact solution for the Ising model in the presence of an external 
magnetic field. Such a solution exists for the one-dimensional analog of the model, although it undergoes no 
phase transition at finite temperature, having very little to contribute to even a qualitative analysis of the 
behavior of real substances. 

Notwithstanding, in the presence of uniform infinite-range interactions phase transitions are induced 
even in one-dimensional systems 0. Furthermore, as pointed out by Hemmer and Stell 0, |8|, for lattice 
and continuous fluids in which the potential has both a long-range attraction and a short-range repulsion, 
plus a hard core, up to two phase transitions are expected. This is indeed verified for the Id lattice gas with 
nearest- neighbor repulsion and long-range attraction R [f0| , described by the Hamiltonian 

N 4:1 N 

H = -4 J n 3 n J+ i - — ^ mn h (1) 
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where J < 0, I > and the occupation numbers iij S {0, 1} represent the absence or presence of a particle 
in the jth site along the chain. Such model is mathematically equivalent to the Ising chain with short- and 
long-range interactions in an external field |Tl| , which in turn is equivalent to a linear chain approximation 
for higher dimensional models. Depending on the ratio of the interaction parameters, the system can present 
one or two phase transitions, and in the latter case a triple point may be present The three possible 

phases have zero-temperature densities given by p — (gas phase), p — \ (liquid phase) and p = 1 (close- 
packed phase). However, the symmetry between unoccupied and occupied cells (hole-particle symmetry) 
|jl2f makes the critical temperatures for both transitions the same, differently from what happens in the 
analogous continuous model (|, [k| . 

On the other hand, for random versions of the model, by looking at ground state properties one can 
see that the presence of disorder in the nearest-neighbor interactions induces both the breaking of the hole- 
particle symmetry and the appearance of additional phase transitions at T = 0. Due to the presence of 
the long-range attraction those new transitions should persist at low but finite temperatures. It is then 
worthwhile examining disordered versions of the model and look for what new interesting properties may 
arise. 

The main purpose of this paper is to study the phase diagram of the Id lattice gas with random annealed 
nearest-neighbor interactions and uniform long-range interactions. This particular disordered version of the 
model [j] has the advantage of being exactly soluble. Although of difficult physical realization, it might be 
useful in studying adsorption problems, if one wishes to consider situations in which the presence of adsorbed 
particles induces diffusion of crystalline atoms in or near imperfect surfaces. 

The outline of this paper is as follows. In order to clarify one subtle aspect of the mathematical treatment 
of annealed models, a preceding discussion is presented. Section [I] discusses some theoretical aspects of 
the statistical mechanics of disordered systems. The d-dimensional lattice gas with an arbitrary uniform 
interaction potential plus random short-range interactions is treated in Sec. following the approach of 
Thorpe and Beeman for the Ising model [Q. The calculation of thermodynamic functions for this class 
of systems is presented, and a general method for calculating the equation of state is proposed. Then the 
one-dimensional lattice gas with random nearest-neigbor interactions and uniform long-range interactions is 
considered in Sec. ^. The main features of the behavior of the model are presented and discussed when the 
nearest-neighbor interactions are selected from a bimodal distribution. In the final section the conclusions 
and main results of the paper are summarized. The equivalence between the model and the random Ising 
model is presented in Appendix A, and in Appendix B the ground state properties of the model are discussed. 



1 Disordered systems. 

Disordered systems are characterized by two kinds of variables: dynamical variables, with mean relaxation 
time t/j, and "structural" variables, related to the randomness of the system and with mean relaxation time 
r c . As introduced by Brout there are two limits in which disorder problems can be formally treated, 
namely, the annealed and the quenched limits. The latter is appropriate to systems where disorder relaxes in 
a much longer time scale than dynamical variables (r c t/j), and so disorder variables can be considered as 
effectively frozen. The opposite situation, in which disorder and dynamical variables fluctuate in the same 
time scale (r c rj th), corresponds to the annealed limit. Since exact solutions can be found for some models, 
in these cases it is worthwhile using the annealed limit as an approximation to a quenched system, as it will 
be further discussed. 
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As shown by Mazo , the free energy of a random quenched system is given by 

F q = -fc B T^p( K )lnQ( K ) + A :B T^p( K )lnp( K ), (2) 
{«} W 

where ks is Boltzmann's constant, T is the absolute temperature of the system, p{n) is the fixed probability 
of occurrence of the random configuration {k}, Q(k) is the partition function of the system for that particular 
configuration, and the summations run over all possible configurations {k}. The second summation in Eq.|| 
is interpreted as due to an entropy of mixing of the random variables. It is in general ignored, since the only 
way to measure an entropy is to change it, which is not possible in this case because the probabilities p(n) 
are rigid by assumption. Irrespective of that term, in general one finds that, for almost all thermodynamic 
functions obtained from the free energy by differentiation, the average value X q of the function can also be 
calculated as an average of the corresponding values X(k) over all the disorder configurations, 

X q = ^ P (k)X(k). (3) 

M 

When lattice models are considered, disorder is usually represented by site or bond random variables. 
For example, in problems of adsorption of gas molecules on a crystal surface one may want to consider 
sites with random potentials or random interactions between particles in different adsorption sites. For 
simplicity, consider a lattice of N sites (or bonds) . Suppose that some random variable Kj can take one of 
M values {-Ki, K2, ■ ■ • , Km} at a given site (or bond) j with probability eft, such that 1i = !• Since 

that probability is independent of the values of the other random variables in the system, the probability of 
a particular configuration {k} = {m, K2, ■ ■ ■ , kjv} is given by 

M 

p(K)=n??*' w 

i=l 

where rii is the number of sites (or bonds) to which the value K j of the random variable is associated in that 
given configuration {k}. Obviously the n^s obey the constraint 

M 



5> = 7V. (5) 



i=l 

It is possible to show that for this case the entropy of mixing in Eq.|| takes the form 

M 

- k B ^ p{n) \xip{k) = -Nk B ^ gi In Qi, (6) 
{«} f=i 

which is proportional to the volume of the system. So, if one keeps the entropy of mixing term and tries to 
calculate (in appropriate units) the pressure of a quenched lattice gas the result is 



P dF « k T 



^2 p{k)P(k) - ^2 H ln Qi 



(7) 
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where P(k) is the pressure calculated for a particular configuration {n}. Clearly Eq.^is a physically incorrect 
expression for the pressure, since in the zero density limit the first term vanishes, while the second term is 
positive. Furthermore, in the non-random limit (Ki = K for all i) all configurations are identical, but the 
entropy of mixing still contributes to the "pressure" . Finally the entropy of mixing diverges in the continuum 
limit (M — > oo, qi — > 0) Jl4| , Thus the correct expression for the pressure must lack the entropy of mixing, 
taking the form 

P = fc B T^>( K )P( K ), (8) 
{«} 

which is the expression used in previous studies of random lattice gases [|l7| [l8|. The need to drop out the 
entropy of mixing in the calculation of the pressure of a quenched fluid was also pointed out by Singh and 
Kovac [fl9[ on the basis of the interpretation of that term as an information entropy, introduced by Sobotta 
and Wagner Q . 

For most systems the quenched limit is more realistic than the annealed limit, since in general t c 3> r^, 
as it happens in spin glasses, for instance. However, the exact mathematical treatment of the quenched 
limit is very difficult, due to the fact that the calculation of the free energy involves averaging logarithms of 
partition functions. On the other hand, the annealed free energy can be written in the form 

F a = -k B T\nQ a (9) 

where Q a denotes the average of the partition function over the disorder configurations. This is often 
calculated with the constraint that the thermal averages of the random variables obey some prescribed 
distribution, which is realized by the introduction of Lagrange multipliers playing the role of pseudo-chemical 
potentials. This implies that in the annealed limit the disorder variables are adjusted so as to minimize the 
free energy, a procedure which singles out a subset of the disorder configurations and introduces correlations 
between the disorder variables themselves. So, the annealed limit is in general not a good approximation 
to the quenched limit. Nevertheless, Morita |2l|] showed that a quenched system can be represented by a 
fictitious (annealed) equilibrium system subject to an additional potential, composed of an infinite series of 
Lagrange multipliers adjusted by the (quenched) n- "particle" correlation functions (n — 1,2,3, . . .) for the 
random variables. Indeed, for various Ising models it can be verified that by controlling correlations between 
pairs of random variables (n — 2) results for annealed systems approach those of the corresponding quenched 
systems |^2|, |23| . This treatment introduces mathematical difficulties which are of course increasingly greater, 
and so this work will analyze the usual case in which one tries to control only the concentration of the random 
variables (n = 1). 

In the next section it will be shown that an entropy of mixing term can also be identified in a large class 
of annealed models, and must likewise be discarded in order to produce meaningful results. 



2 Lattice gas with random short-range interactions: annealed 
limit. 

Consider a iV-site d-dimensional lattice gas model whose particles interact via an arbitrary uniform (non- 
random) potential plus short-range interactions satisfying the annealed probability distribution 

M 

p( K v) = 5^ - J k ), (io) 
fe=i 
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for all pairs of nearest-neighbor sites i and j. The Hamiltonian of the model can be written as 

H = H^+Y,H ih (11) 



<hj> 



where H^ u ' is the uniform term, the summation runs over all nearest-neighbor pairs of sites and 



M 



k=l 



(12) 



For any site i, variable rii = if the site is empty or n,; = 1 if the site is occupied by a particle. 

Following the approach of Thorpe and Beeman for the Ising model [Q, this situation can be formally 
treated by writing K%j in the form 

M 

K ij = ^ ' tjj kJk, (13) 
k=l 

and introducing M pseudo-chemical potentials to control the averages of the random numbers which 
for all are subject to the constraints 



M 



ty.fc S {0, 1}, }J ijt k = 1, and (*y,fc) = 



fc=i 



The annealed grand partition function of the system is given by 



Z = ^2' exp 

it} 



M 



k=l <i,j> 



J^exp 

{«} 



N 



3=1 



(14) 



(15) 



where the primed summation runs over those configurations satisfying the first two constraints in [w]. After 
performing the partial trace over the random variables one obtains 



Z = 2_, ex P 



A' 



j=i 



<i,j> lfe=l 



Due to the fact that rii € {0, 1} for all sites the term in braces in Eq.[l6] can be expressed as 



M 



J2 exp [P (& + JfcTunj)] = Ae Kl 



k=l 

if A and K are defined by the expressions 

M 

A = w k and e K = wy 

k=l k=l 

with = exp(/3^). Equation [l6]can now be written as 

N 



M I M 

0J k 



k=l 



Z = A N q /2j2exp 

{n} 



3=1 <i,j> 



= A Nq / 2 Q(K), 



(16) 



(17) 



(18) 



(19) 
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where q is the coordination number of the lattice and Q(K) is just the grand partition function of a regular 
lattice gas with short-range interaction energy K//3 plus the uniform interaction potential. 

The fugacities Wk are eliminated by imposing for all nearest-neighbor pairs (i,j) and all values of k the 
condition 

. . Wk dlnZ 

{Uj,k) = Qk = Qk, (20) 



which can be rewritten as 



where 



din A 

dw k 



e(K)w k 



N dw k 
dK 



qk, 



e(K) = {UiUj) 



31 K 



dwk 

2 dlnQ(K) 
Nq dK 



(21) 



(22) 



is the nearest neighbor pair correlation function for the regular lattice gas with interaction Kj (3. By rear- 



ranging terms in Eq.|2l| and imposing the condition Y2k=i qk — ^ °ne obtains 

M 



E 



Qk 



I coth [\{K- pJ k )] + 1 - 2e(K) 



0, 



(23) 



which determines K if one is able to calculate e(K). Notice that in general the effective interaction K/(3 
depends on the temperature. 

The density of the gas is set by the condition 



P = 



d\nZ _ z dlnQ(K) 



N dz N dz 

and from this equation it is possible to eliminate the fugacity z. 
The annealed free energy is given by |24 



(24) 



fa = -k B T 



1 M 

— In Z - ^ qk In w k - p In z 



k=l 



(25) 



which after straightforward calculations can be written in the form 

C M 



M 



fa = fx -k B T\J2 Ik In [1 - e(K) (l - eP J "- K )] - ^ q k hxq k 



k=l 



where 



fx = -k B T 



— In COR") -pin* 



(26) 



(27) 



is the free energy of the regular system. The last summation on the right hand side of Eq.|2^ is an entropy 
of mixing term, identical to that in Eq.^, and must be dropped as discussed in Secjl]. With that correction, 
the free energy per lattice site takes the form 



M 



f = f K - k B Tj2qk In [1 - e(K) (l - e^ K )] , 



(28) 



fc=i 
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and the pressure of the system can then be calculated from the thermodynamic relation 



P = PH-f = k B Tpinz-f=> (29) 

=> P = k B T 1 1 In Q(K) + » In [l - e(K) (l - e^~ K )] j , (30) 

the first term being just the pressure of the regular lattice gas with effective interaction K/j3. In the non- 
random limit (Jj. = J for all k) Eq.23 gives K — (3 J and according to Eq.^the pressure is given by 



P = P K =^UQ{K), (31) 

which is obviously the correct expression. 

From the results of this section it is then clear that, as in the quenched limit, the entropy of mixing term 
can be made explicit in the free energy of any lattice gas with annealed short-range interactions, and that 
the correct equation of state can be obtained by discarding that term. This is a fairly general result and will 
be applied in the next section to the Id lattice gas with random nearest-neighbor and uniform long-range 
interactions. 



3 One-dimensional lattice gas with random short-range interac- 
tions and uniform long-range interactions: annealed limit. 

The existence of exact solutions [l(J for the uniform model [j] allows one to exactly solve the model in the 
presence of random annealed nearest-neighbor interactions, as shown in the previous section. 

The Hamiltonian of a Id lattice gas with uniform long-range interactions and random nearest-neighbor 
interactions is given by 



/v 



H 



3 = 1 



41 ^ 
— > 



where the nearest-neighbor interactions Kj satisfy the annealed distribution 



p( K i) = ^QkS(K>j - Jk, 



(32) 



(33) 



fc=i 



The thermodynamic properties of this random model can be obtained from those of a reference system 
without the long-range interactions. In particular the Hclmholtz free energy is given by B 



f(p,T) = CE f{p,T)-U? 



(34) 



where /(p, T) is the free energy of the reference system as a function of the density p and temperature T, 
while "CE" denotes the convex envelope of the function in square brackets. Correspondingly, the pressure 
P and chemical potential p are given by fl 



P = MC 



P - 4 Iff 



p = fx - 81 p, 



(35) 
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where "MC" denotes the Maxwell construct of the function in square brackets and 

df 

P = PP>~ f\P, T) and A = g- 



(36) 



are the corresponding functions for the reference system. 

From the results of the previous section, Eqs.|2^ and the free energy / and pressure P of the reference 
system without long-range interaction are given in terms of the corresponding functions fx and Pk for a 
regular model with effective nearest-neighbor interaction K/ (3. The grand-partition function of this regular 
model is given by 

N N 

0ji rij +H'Y^ rij n j+1 1 , (37) 



Q(X) = ^exp 
{«} 



and by using the transfer matrix method with cyclic boundary conditions one easily finds 



I3P K = In Q{K) = In 1 <j 1 i< 



z„4K 



fl ze iK ) 



45 



1/2 



(38) 



where 5 = exp (/3/i) . The density of the system is then given by the relation 



OPk 



[(1 



~ze AK ) e 4K 



2] 



■45 



-1/2 



1 + Se 4K 



(l~Se 4K ) +42 



1/2 



(39) 



from where the fugacity z of the reference system can be determined, and the nearest-neighbor pair correlation 
function of the regular system is given by 



e(K) 



d{(3P K 
~dK~ 



ze 4K 



l-(l~Se 4K ) (l-ze iK )+4z 



-1/2 



1 + ze 



(1 - ze iK ) +45 



1/2 



(40) 



Finally, the pressure of the random system with long-range interactions described by the Hamiltonian in 
Eq.^| is 

C M ~\ 



P = MC < Pk + k B T Qk In [l - e(K) (l 



40J k -4K 



)]-4I, 2 



(41) 



k=l 



with K determined from the solution of Eq.|23|. So, the complete equation of state can be obtained from 
Eqs.|8H4l]. 

An alternative approach can be used to solve this problem, namely, the mapping onto the random Ising 
model. This is discussed in Appendix [| and has the advantage of making evident the breaking of hole-particle 
symmetry. Furthermore, it provides a way of avoiding the use of the Maxwell construct by previously locating 
free energy minima, which directly indicate phase transitions. 

The results above are applied to the case in which the nearest-neighbor interactions are selected from 
the bimodal distribution 

p(Kj) = pS(kj - J A ) + (l-p)S(Kj - J B ), (42) 
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where the interactions Jb are assumed repulsive (Jb < 0). This is done in order to introduce competition 
effects, which, as in the pure model, are responsible for complex behavior. In the absence of competition 
between short- and long-range interactions the system presents at most one phase transition. The ground 
state properties of the model, discussed in detail in Appendix ||, indicate that at T = and for J a ^ 
the system undergoes at most three phase transitions. For J a > 0, when there is competition between 
the nearest-neighbor interactions, at fixed pressure the system can present itself in structures characterized 
by p — (gas phase), p = p, p — ^(1 + p) (liquid phases), and p = 1 (close-packed phase). For J a < 
all nearest-neighbor interactions are repulsive, and the stable phases are those for which p = (gas), 
p = |, p = |(1 + p) (liquid phases) and p = 1 (close-packed). For J a = the system undergoes at most 
two transitions, associated with those of the pure model |§, [h]], the liquid-phase density being shifted to 
p=\{l+p). 

The characterization of the phases p = i, p = p and p = i(l + p) as liquid phases is based on the 
instability of the corresponding structures. As in the uniform case, the p — \ phase at T — is stable at 
certain pressures because particles are forced by nearest-neighbor repulsion to occupy every other site of 
the lattice. At any finite temperature, thermal fluctuations overcome the effect of the repulsion, inducing 
local fluctuations of the density and destroying the structure. The stability of the p = p and p = ^(1 + p) 
phases at T = 0, on the other hand, is related to the existence of interaction domains (Appendix ||). These 
are a consequence of the interactions mobility in the annealed limit and are formed as particles are added 
to the system at fixed volume in order to minimize the free energy. In the p — p phase all sites in the J a 
domains are occupied, while all those in Jb domains are empty. For p = i(l +p) all sites in J a domains and 
every other site in Jb domains are occupied. However, interaction domains are also destroyed by thermal 
flucuations at any finite temperature, making the corresponding structures unstable. 

In the following subsections the properties of the model at all temperatures are discussed for the various 
possible cases. The relevant parameters are renormalized by \Jb\, 

,_Ja _ I „_k B T p . = P (43) 



| Jul' \JbY ~ \JbV \Jb[ 

In passing it should be mentioned that, for repulsive long-range interaction (/ < 0), the system does not 
undergo a phase transition irrespective of the sign of the short-range interaction. 



3.1 Case J A = 0. 

For dilute short-range interactions there are at most two phase transitions, whose critical temperatures 
are different for p £ (0,1) because of the breakdown of hole-particle symmetry, as already stated. The 
critical temperature T c of the gas-liquid (G-L) transition is always higher than T' c , the critical temperature 
of the liquid-close-packed (L-CP) transition. This is in agreement to what is verified for the corresponding 
non-random continuum model [13| , in which there is no hole-particle symmetry. 

For p ^ 1 there is always a range of values of a in which one or two triple points exist, as can be seen from 
the reentrant behavior of the triple point temperature (T t ) curves in the T x a diagrams of Figs. l(a)-(c), 
for p = 0.1, p = 0.25 and p = 0.5, respectively. The critical temperature of the G-L transition coincides 
with the triple point temperature for a — a ctl whose dependence on pis shown in Fig. 1(d). It can be seen 
that a c t approaches at r — 3.1532 as p — > 0, while for p — » 1, where there is no short-range repulsion, a c t 
approaches ao = 2. This is the value of a for which T t vanishes for all p when J a — 0. 

Examples of the behavior of the isotherms and phase diagrams for p = 0.1 are shown in Figs. 2-4. In 
Fig. 2, for a = 1.5, there are two-phase transitions even at T = and there is no triple point. For low 
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temperatures the derivative of the L-CP transition pressure with respect to T is negative. This is related to 
the destruction by thermal fluctuations of the interaction domains built at T = 0, which makes the effective 
short-range repulsion weaker, lowering the pressure needed to pack the particles. In Fig. 3, for a = 2.2, there 
is only one phase transition (G-CP) at T — 0, although two transitions occur for a range of temperatures 
above the triple point. In Fig. 4, for a = a ct ~ 2.4297, the critical temperature of the L-CP transition 
coincides with the triple point temperature, and for a > a ct there is only one phase transition. The behavior 
of the system when there are two triple points will be exemplified in the next subsection. 

3.2 Case J A > 0. 

When there is competition between short-range interactions even a small amount of randomness is responsible 
for interesting behavior, as can be seen in the p x T projection of the coexistence surfaces shown in Figs. 
5(a) and (b) for J a = — Jb {S = 1), p = 0.07 and two slightly different values of a. In both cases there 
are three phase transitions at T = 0, defining in decreasing density the close-packed, the two liquid (Li and 
L 2 ) and the gas phases. The lowest critical temperature corresponds to a L 2 -G transition. Above that there 
is an intermediate range of temperatures in which the system presents only two transitions. For a = 1.13 
that range is limited by a triple point, immediately above whose temperature there is again three phase 
transitions, as shown in Fig. 5(a). But for a — 1.129 that triple point is replaced by another critical point, 
giving rise to a coexistence "bubble" , as can be seen in Figs. 5(b). The T x a diagram for 6 = 1 and p = 0.07 
is shown in Fig. 6. When the long-range attraction is sufficiently weak the system can present three phase 
transitions, whose critical temperatures are shown as solid lines in the figure. Triple points (broken lines) 
are present for intermediate values of a. The breaking and the reentrant behavior of the T c curve, evident in 
the inset of Fig. 6, is related to the appearance of coexistence "bubbles". That unusual behavior is observed 
only for values of p between p ~ 0.06 and p ~ 0.08. 

The system may also present a quadruple point, i.e. coexistence of four phases, if J a > — Jb {S > 1). 
This is exemplified for 5 — 2 and p — 0.1 in the T x a diagram in Fig. 7, where four triple point curves 
(dotted lines) intersect at the quadruple point P q . Figures 8 and 9 show P xT and pxT diagrams for 5 = 2, 
p = 0.08 and values of a around a q , illustrating the possible existence of two triple points or a quadruple 
point. 

3.3 Case J B < J A < 0. 

As detailed in Appendix B, when all short-range interactions are repulsive the possible structures at T = 
correspond to p = 1 (close-packed phase), p = |(1 + p) (liquid- 1 phase), p = | (liquid-2 phase) and p = 
(gas phase). 

When \Jb\ is much larger than \Ja\ the behavior of the system is expected to be similar to the dilute 
case {J a = 0). This is confirmed for the case J a = \\ Jb {5 = — jq), as shown in the 9 x a diagrams in Figs. 
10(a)-(c) and the a x p diagram in Fig. 10(d). Except for the presence of three phase transitions for small 
values of a, the behavior is similar to that observed in Fig. 1. Isotherms and phase diagrams are shown 
in Figs. 11 and 12 for p — 0.5, a = 0.8 and p = 0.75, a = 0.8, respectively. As expected for S > — \ (see 
Appendix B) the two-transition regime at T — consist of the p = 0, p = ^(1 + p) and p = 1 phases, 
as in the dilute case. For small values of the long-range interactions, the great energy difference between 
short-range interactions leads to the high pressure needed to close-pack the system, when it is necessary 
to occupy nearest neighbor sites with Jb bonds. At low pressure the occupation of those sites is virtually 
forbidden, giving rise to the behavior shown in Figs. 11 and 12. As the value of a is increased, the system 
becomes more similar to the dilute case, as can be seen in Figs. 10(a)-(c). 
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When \ \Jb\ < \Ja\ < %\Jb\ {— § < <5 < —5) both two-transition regimes at T = 0, described in 
Appendix B are possible, with the liquid-phase corresponding to p = A or p = i(l + p) if p is less or greater 
than p' = (1 + 2<5)/(l + 5), respectively. At finite temperatures the case p — p' contains ingredients of both 
regimes. The 9 x a diagram for the case J a — {S — — § ) and p = p' = | is shown in Fig. 13. Notice 
that both triple point lines touch the a axis at the same point, indicating the absence of a two-transition 
regime at T = 0. Figure 14 shows isotherms and phase diagrams for the same case and a = 1.45, where a 
triple point is observed, corresponding to the coexistence of the close-packed and the two liquid phases. The 
L1-L2 transition line has a positive temperature derivative, as shown in Fig. 14(b), differently from what is 
observed for 5 > 0, as shown in Fig. 8(a). 

For p < p' and in the two-transition regime, the T = liquid-close-packed transition pressure is given by 
P = —4[pJa + (1 — p) Jb + 1], the arithmetic mean of the corresponding results for the pure cases p = and 
p=l. For p > p' the system can present a quadruple point, as observed for the 5 > 1 case. 



4 Conclusions. 



The Helmholtz free energy of any lattice gas model with random short-range interactions in the annealed 
limit was shown to present an entropy of mixing term. As shown for the quenched limit, that term has to 
be dropped from the free energy in order to give correct physical expressions for the pressure. 

As an application, the one-dimensional lattice gas with uniform infinite range interactions and random 
short-range interactions satisfying bimodal distributions was exactly solved. Explicit results were presented 
for the cases J a = 0, Ja = —Jb > 0, Ja = —2Jb > 0, J a = JqJb < and J a — \ Jb < 0, where there is 
competition between short- and long-range interactions. In general, the system can present at most three 
phase transitions at a fixed temperature. The phases were identified as close-packed, gas and two liquid 
phases, the latter being characterized at low temperatures by the existence of interaction domains. A very 
important effect of disorder was the breaking of the hole-particle symmetry of the pure model, which induces 
different critical temperatures for the transitions. Under certain conditions phase diagrams present negative 
thermal derivatives of some transition pressures, a phenomenon observed in nature in water, for instance. 
The existence of triple and quadruple points was demonstrated, and for Ja > and small values of p the 
existence of coexistence "bubbles" was observed. 



Appendix A. Equivalence between the random lattice gas and the random 

Ising model. 

Here an alternative treatment for the problem of the random lattice gas discussed in Sec. || is presented. 
By using the transformation Q 

n j = (1-^0/2, <r s = ±l, (44) 
it can be shown that the annealed grand partition function of the random lattice gas, 



Z = ex P 



M N 

fe=i j=i 



E ex p 



N 



J'=l 



(45) 



can be written as 



Z = exp 



k=l 



(46) 
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where / = (51, h = (3h — —21 — \(3^l and 



Z„ = E exp 

{t,a} 



N M j N N 

/? E E l 3- k + J k<7j<rj+i) + ]y E ai(7 i + E ^. 



3=1 k=l %,3 = 1 

is the annealed grand partition function of a random Ising model with 

ik = S,k + Jk 

and 

M 

hj =h + 20j2^k-tj,k] J k' 



3 = 1 



(47) 



(48) 



(49) 



k=l 



Notice that the distributions of the random fields hj and the short-interactions Kj are correlated. The 
presence of the random terms in the effective field hj acting on the spin in the jth site breaks the h — > —h, 
{a} — > {— a} symmetry present in the uniform Ising model (for which J k = J for all k). The effect of this 
on the behavior of the gas is that, unlike the pure model, the hole-particle symmetry is also broken in the 
presence of any randomness, and so all transitions have different critical temperatures. This is in contrast 
to the lattice gas with particles interacting via softened core uniform potentials {?], |§| , where multiple phase 
transitions are also present, but whose critical temperatures are not all different because the symmetry is 
still present. 

The correlation between the random fields and random interactions allows one to perform the partial 
trace over the disorder variables, mapping the system onto a regular Ising model with an effective interaction 
KJ (3 and a interaction dependent effective field. Rewriting Eq.|^] as 



n r m 

{a} 3 = 1 U=l 



PJk (OjOj+\ - (Tj - (Tj+l) + f3£k 



where 



Z(a) = exp 



j N / M \ N 

— e a i a 3 + [ ~ h + 2 @ E qk Jk ) E a i 



i,3=l 



k=l 



3=1 



and performing a decimation 

N ( M 

n e«p 

3=1 U=l 

with 



(3J k {UjUj+i - Uj - Gj+l) + I3£, k 



A* 



t M 

E 

\k=l 
/ M 

E 

\k=l 



w k e 



3/3 J k 



' M 

E 

\k=l 



w k e 



-0J k 



3p.J k 



t M 

E 

\k=l 



w k e 



(50) 
(51) 

(52) 

(53) 
(54) 
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and 

w k = exp (Pik) , 
one can express the annealed grand-partition function Z a as 

N ? N 



Z a = A N J2 ex P 



N j N N 



3 = 1 



i=i 



where 



h = h + 2 ^Io + pY^qkJu -K^j 



(55) 



(56) 



By eliminating A and the fugacities Wk the annealed free energy of the random Ising model (as a function 
of the magnetization) can be written as 



M 



fa 



fc=l 



qk In qu 

M 

qk hi 



k=l 



[3 - e(K)]e^ K -' 3J ^ + [1 + e(K)}e- 2 ( K -P J ^ ' 



in which once more there is an entropy of mixing term and 



f(K,a) 



-knThi < e 



cosh h + ^sinh 2 h + e 



-AK 



1/2' 



1(7 



(57) 



(58) 



is the free energy of an Ising model with average magnetization a in a uniform field and uniform short- 
(K/P) and long-range (J) interactions flO|| , K being determined from the solution of the equation 



M 



^ 2coth [2 (# - 13 J k )] + [1 - e ff (iT)] 
The nearest-neighbor spin-spin correlation function is explicitly given by 



0. 



(59) 



e a {K) = (aja j+1 ) K = 1 



2e- iK (sinh 2 h + e^ K 
cosh h + (sinh 2 h + e~ iK 



-1/2 



1/2- 



(60) 



By inverting the transformation 44 the annealed free energy of the lattice gas can be written as 

M / M \ 



k B T 
N 



+ q k £,k + PP = [h + ^qkJk + I ) + fa 



PPi 



(61) 



fc=i \ fc=i / 

from where the pressure is calculated by explicitly excluding the entropy of mixing term present in f a , 

M 

P = pfi — f g + k B T Qfe In gfe => 
fe=i 
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(M \ M 

ft + J^^Jk + I - + k B T^q k Ing fc> (62) 
fe=i / fe=i 

while the density is given by 

P=^. (63) 

By studying the properties of the Ising model described by the free energy f„ in Eq.[57] the behavior of 
the lattice gas can be determined. The equilibrium state of the magnetic system at fixed temperature and 
external field h is that which minimizes f a . By locating those minima and using Eqs.|62| and [33] the equation 
of state of the gas is exactly determined. This procedure automatically removes the unstable portions of the 
isotherms and is equivalent to the Maxwell construct. 

Appendix B. Ground state properties. 

At zero temperature the Helmholtz free energy F of a system is identical to its internal energy, as it is 
evident from the thermodynamic relation 

F = E-TS, (64) 

where T is the temperature and S is the entropy. 

The ground state behavior of the system can thus be determined by writing the internal energy E for 
the various ranges of the density p and calculating the pressure using the expression 



P\t=o = MC 



p—-E 
dp 



(65) 



where "MC" denotes the Maxwell construct || of the function in square brackets, which substitutes horizontal 
parts (corresponding to phase transitions) for the unstable portions of the isotherms. 

A fundamental feature of annealed models is that impurities are free to rearrange themselves so as to 
minimize the free energy. For a lattice gas with random interactions this means that, as particles are added 
to the system, interactions of the same kind have a tendency to aggregate, forming domains. From this 
analysis it is possible to find those structures which minimize the internal energy and determine the ground 
state behavior of the model discussed in Sec. ||, which is summarized below. In the following discussion the 
volume of the system is regarded as fixed. 

Case Je < < J a- In this situation, because of the presence of the attractive interactions J a, interaction 
domains are formed even at infinitesimal densities. For sufficiently weak long-range attraction energy J, the 
density can be increased from p = to p = p at zero pressure by adding particles to the J a domains. At 
densities between p = p and p = h(l+p) new particles can be allocated with minimum energy cost at every 
other site of the Jb domains. The pressure is nevertheless raised, because of the smaller volume available. 
Finally, for p > i(l + p) particles can be added only by allocating them in empty sites between two other 
particles in the Jb domain, raising both internal energy and pressure. For strong long-range attraction some 
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of the transitions are suppressed, the various possibilities for the T = isotherms being as follows. 



3 transitions: 

< p<p, P = = Pi 

p<p<±(l+p), P = ApJ A -2p{l+p)I = P 2 

|(1 + p) < p < 1, P = A[p.J A - (1 - p) J B ] - 2(1 + p)I = P 3 

condition: J < min j-^ Js, Ja| 
2 transitions: 

r o<p<i(i+ P ), p = pj 
\ |(i + P ) <p< i, p = p 3 



conditions: 



T ^J A </<2[ T ^J^ 



2 transitions: 

f 0</9< P , P = Pl 

[ p < p < 1, P - 4p(J A -J B -I)=P! 2 

!Ja> —Jb and 
p < Pd and 
"T^ J s < I < Ja-Jb 

1 transition: 

{0<p<l, P = Pi 
p < _p d and I > J A - J B 
or 



conditions: 



p > Pd and / > 2 



Case Jb < Ja < 0. This is the case where all nearest-neighbor interactions are repulsive. For sufficiently 
weak long-range interactions / there is no aggregation of particles at low densities. On the contrary, particles 
can be added to the system at minimum cost by allocating them such that they have no nearest neighbors. In 
general this can be done at zero pressure only if p < |. For greater densities interaction domains are formed, 
and new particles are allocated, by raising the pressure, between occupied sites in the Ja domains, which 
are filled up for p = p+ |(1 — p) = |(1 +p). Adding further particles requires even higher pressures, since 
they must occupy available sites in the Jb domains. Again the existence of the various possible transitions 
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depends on the strength of the long-range interactions. Zero temperature isotherms are described below. 

3 transitions: 

0<p<i, P = = Pl 

I <p< P = -4J A -(l+p)I = P 2 

|(1 + p) < p < 1, P = 4[pJ A - (1 - p)J B ] - 2(1 + p)J = P 3 

condition: J < min j — J^, 4(J^ — Jb) j 
2 transitions: 

r o<p<i(i+ P ), p = A 
\ i(i+ P )<p<i, p = p 3 



conditions: < 



Ja>\Jb and 



^ 2J A -J B 
P ^ Ja-Jb 



p s and 



^J^</<2 



Jb 




conditions: 



2 transitions: 

P = Pi 

P = -A[pJ A + (1 - p) Jb + 1] = P 2 ' 
Ja<\Jb and 
p < p s and 

4( Ja - Jb) < / < -2[p J4 + (1 - p)Jb] 



conditions: 



1 transition: 

{0<p<l, P = Pi 
p < p s and / > —2\pJ^ 
or 



p > p s and I > 2 



l+p Ja 



(i-p)Jb] 

- Jb 
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Figure captions 

Figure 1. (a-c) 9 x a diagram, where 9 = ksT/ |Js| is the renormalized temperature and a = I / \ Jb\, 
for dilute short-range interactions ( J a = 0) with probabities p = 0.1, p = 0.25 and p = 0.5, respectively. 
Continuous curves represent critical temperatures (0 C and 0£) 01 the phase transitions, while triple point 
temperatures (9 t ) are shown as dotted curves. Curves for 9' c and 9 t coincide at point P ct , whose a coordinate 
(ad) is shown in (d) as a function of p. 

Figure 2. (a) Isotherms P* x p (continuous curves), where P* = P/\Jb\ is the renormalized pressure, 
and coexistence curves (dotted curves) for dilute short-range interactions (Ja = 0), p = 0.1 and a = 1.5; 
(b) Phase diagram P* x 9 showing close-packed (CP), liquid (L) and gas (G) phases. Continuous curves 
correspond to first-order transitions terminating at the critical temperatures 9 C and 8' c . 

Figure 3. (a) The same as in Fig. 2(a) for a = 2.2; (b) The same as in Fig. 2(b), but with a triple point 
indicating coexistence of the three phases. 

Figure 4. (a) The same as in Fig. 2(a) for a ct ~ 2.4297; (b) Phase diagram showing the coincidence of the 
triple point (temperature t ) and the critical point of the liquid-closc-packed (temperature 8' c ). 

Figure 5. p x 6 projections of the coexistence surface for S = 1 (6 = Ja/ \ Jb\), P = 0.07 and (a) a = 1.13 
and (b) a — 1.129. The triple point present in (a), whose temperature is 9 tl is suppressed in (b), where a 
coexistence "bubble" can be observed, as indicated by the arrows. 

Figure 6. 9 x a diagram for S = 1 and p = 0.07. The inset shows the breaking and reentrant behavior 
of the 9 C curve, as well as the existence of the new points P" t and P" t ', extremes of a new triple point line 
(dotted). The 9' c curve is ommitcd from the inset for clarity. 

Figure 7. 9 x a diagram for 5 = 2 and p = 0.1. The four (dotted) triple point lines intersect at the 
quadruple point P q . 

Figure 8. Phase diagram for 5 = 2, p = 0.08 and (a) a = 2.3, (b) a q ~ 2.3277 and (c) a — 2.35. Two triple 
points are evident in (a) and (c), while the quadruple point is present in (b). In all cases the axes intervals 
are < 9 < 2.4 and < P* < 0.8. 

Figure 9. px 9 (9 = k B T/ \ Jb\) projections of the coexistence surfaces for S = 2, p = 0.08 and (a) a = 2.3, 
(b) a q ~ 2.3277 and (c) a — 2.35. The broken lines indicate the triple point temperatures in (a) and (c), 
and the quadruple point temperature in (b). 

Figure 10. (a-c) 9 x a diagram, where 9 = k B T/ \ Jb\ is the renormalized temperature and a = 1/ \Jb\, for 
the case Ja — Jq Jb < (S = — j^) with probabities p = 0.25, p = 0.25 and p = 0.5, respectively. Continuous 
curves represent critical temperatures (9 C and 9' c ) of the phase transitions, while triple point temperatures 
(9 t ) are shown as dotted curves. Coordinates a ct and a' ct of points P ct and P' ctl as well as values ot2 and a 3 
at which triple point temperatures become zero, are shown in (d) as a function of p. Notice that, while a ct 
goes to atr — 3.1532 as p — * 0, a' ct goes to \S\ at r — 0.9459 as p — > 1, as expected. 

Figure 11. (a) P* x p isotherms (solid lines) and coexistence curves (dotted lines) for 5 = —0.3, p = 0.5 
and a — 0.8; (b) P* x 9 phase diagrams showing the various transition lines and critical temperatures; (c) 
p x 9 projections of the coexistence curves shown in (a). 

Figure 12. The same as in Fig. 11 for p = 0.75. Notice in (b) and (c) the existence of a triple point. 

Figure 13. 9 x a diagram for 5 = —0.6 and p = p^ = 0.5, showing the critical temperatures (solid lines) 
and triple point temperatures (dotted lines). 
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Figure 14. (a) P* x p isotherms (solid lines) and coexistence curves (dotted lines) for S = —0.6, p = 0.5 
and a — 1.45; (b) P* x 9 phase diagrams showing the various transition lines, critical temperatures and the 
triple point. 
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Figure 1. A. P. Vieira and L. L. Gongalves. 
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Figure 2. A. P. Vieira and L. L. Gongalves. 
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Figure 3. A. P. Vieira and L. L. Gongalves. 
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Figure 4. A. P. Vieira and L. L. Gongalves. 
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Figure 5. A. P. Vieira and L. L. Gongalves. 
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Figure 6. A. P. Vieira and L. L. Gongalves. 




6 



Figure 7. A. P. Vieira and L. L. Gongalves. 
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Figure 8. A. P. Vieira and L. L. Gongalves. 




8 



Figure 9. A. P. Vieira and L. L. Gongalves. 
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Figure 10. A. P. Vieira and L. L. Gongalves. 
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Figure 11. A. P. Vieira and L. L. Gongalves. 
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Figure 12. A. P. Vieira and L. L. Gongalves. 
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Figure 13. A. P. Vieira and L. L. Gongalves. 
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Figure 14. A. P. Vieira and L. L. Gongalves. 
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